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We use a dispersion relation in conjunction with the operator product expansion (OPE) to derive 
model independent sum rules for the dynamic structure functions of systems with large scattering 
lengths. We present an explicit sum rule for the structure functions that control the density and spin 
response of the many-body ground state. Our methods are general, and apply to either fermions 
or bosons which interact through two-body contact interactions with large scattering lengths. By 
employing a Borel transform of the OPE, the relevant integrals are weighted towards infrared fre- 
quencies, thus allowing for greater overlap with low energy data. Similar sum rules can be derived 
for other response functions. The sum rules can be used to extract the contact parameter introduced 
by Tan, including universality violating corrections at finite scattering lengths. 
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I. INTRODUCTION 



The theory of non-relativistic spin-1/2 fermions near the unitarity hmit (with 5- wave scattering length a oo), 
described by the Lagrangian 

plays a role in several areas of physics, ranging from nuclear physics to atomic and condensed matter physics. Due to 
the recent experimental progress in manipulating gases of trapped cold atoms tuned to a Feshbach resonance, it has 
become a problem of relevance to obtain first principles predictions for the many-body collective behavior of fermionic 
systems near unitarity. 

Although the dynamics in the limit a — >■ oo is universal, the lack of the textbook small expansion parameter, na^, 
with n the atomic number density, makes it difficult to derive model independent predictions for the interacting 
quantum many-body ground state of the system. However, starting with recent work of Tan [H-Q, it has become 
apparent that many observables of the system are governed, in the a — > c» limit, by a set of universal relations which 
hold true in any quantum state of Eq. ([T|), and which can be fixed by studying the short distance/time scale behavior 
of the system. These observables are summarized in the recent review H, and besides those found in include rf 
spectroscopy [1,0, photoassociation [3, 13], static [9j] and dynamic [lol. structure factors, and hydrodynamics [T]|. 

From the point of view of quantum field theory, the existence of universal relations which hold in any quantum 
state are a consequence of of an enhanced non-relativistic conformal (Schrodinger [l2| - [l5j ) symmetry in the ultraviolet 
(corresponding to the limit a — > oo), together with the operator product expansion (OPE). The OPE was originally 
introduced by Wilson 16] in the context of strong interaction physics, and the connection between the OPE in non- 
relativistic quantum field theory and universal relations such as those of [H-Q was first made explicit in the work of 
Braaten and Platter which reproduced and extended the original Tan relations using the OPE. More recently, 

refs. ^,1^ showed that the OPE predicts that the asymptotic large frequency behavior of dynamic structure functions 
is also controlled by the Tan contact parameter C. 

By itself, the OPE only controls the properties of short distance/time observables, e.g., tails of distributions in 
the large frequency or momentum limit. However, when used in conjunction with dispersion relations it is possible 
to obtain indirect information about low energy physics, in particular the spectrum of low energy excitations about 
the ground state, from the OPE. In this paper, we will use the OPE in conjunction with a set of weighted ("Borel 
transformed" ) dispersion relations to derive sum rules that probe the infrared properties of many particle systems 
whose dynamics is described by Eq. ((T]). The sum rules relate weighted frequency integrals of retarded correlators to 
thermodynamic properties of the many-body ground state, in particular the number density, energy density, and the 
Tan contact C. These weighted integrals depend on a free parameter ojq which, if taken small enough, are dominated 
by low energy data, and are therefore sensitive to the spectrum of excitations about the ground state. 

The methods we employ to obtain these sum rules are based on analogous techniques used first used in a relativistic 
setting to compute QCD sum rules. [2^. Our sum rules are applicable both to fermionic systems and to bosons near 
unitarity (e.g., He-4) and include correction to the infinite scattering length limit. The techniques can be applied to 
any response function, but for the sake of illustration we will here consider only the density-density structure factor 
S{uj,q) . 

In secHH we review how the OPE controls the asymptotic behavior of Green's functions, and use this to obtain sum 
rules. We apply this methodology to the structure function S{uj,q) in sec. IIIII The main results are given in Eq. (j70p 
and Eq. ([75]) , and are valid for arbitrary values of the scattering length, but neglect subdominant effects in the OPE 
from the range tq of the two-body potential or from local operators that mediate three-body physics and higher. 
Taking the a — oo limit of our calculation, the results of this section agree with those of fl^. In sec. IIVI we consider 
dynamic spin-density correlators, and find, as expected, that the leading a; — )■ oo asymptotics is entirely fixed in terms 
of the Tan contact C. The sum rules for the spin structure function derived in sec. |IV]are of phenomenological interest 
since they allow for direct access to the contact C from low energy data. Results for spinless bosons are qualitatively 
similar to those obtained for fermion in sec. IIIII and are given in sec. [V] Conclusions and directions for future work 
are discussed in sec. |Vl] 
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II. SUM RULES FROM THE OPE 

A. The OPE and its range of validity 

In quantum field theory, the OPE is the statement that an operator product A{x)B{0), in the limit of short 
distance/time scales, is equivalent to an expansion over local operators at x — 0, 

A(.x)B(0) ^^a(a;)O„(0). (2) 

a 

This relation implied by this equation is strong in the sense that it holds valid inside all matrix elements, i.e. it is a 
statement about the operator algebra itself. Inserted inside a matrix element or correlation function, the OPE relation 
is an asymptotic expansion in the limit where the distance scale x is much shorter than the scales associated with 
external states or with additional operator insertions. 

In a non-relativistic field theory (with d — 3 spatial dimensions), the OPE takes the form 

A{x)BiO) ^ 1^1^°"^"'^"/- (^) ^^-(0), (3) 

In this equation, Aa^b.cx are the scaling dimensions of the operators A{x)^ B{x),0{x) respectively, defined by 

0{Xx,X'^t) = X-'^0{x,t), (4) 

and the coefficient function fa{\x\'^/t) is a calculable function, as we will discuss in detail in sec. IIIII Taking ground 
state matrix elements on both sides of Eq. ^ (or ensemble averages at finite temperature T), one finds that the 
asymptotic short distance or time properties of two-point Green's functions, for instance the response functions of the 
theory, are dominated by the condensates (Oa) of the first few operators of lowest dimension A^. 

In this paper, we will use Eq. ([3]) to determine the asymptotic behavior of response functions and to derive sum 
rules for frequency dependent transport coefficients. Ideally, one would like use the OPE to constrain the retarded 
Green's functions, 

iGniq)^ [ d^xe'^--e{t){[Oix),Om), (5) 



where q'' = (cj, g), and (• • • ) is either a vacuum expectation value (at T = 0) or a thermodynamic average (• • • ) = 
Tr[e~'^'^~^^' • • ■]/Z. In this equation the operator 0{x) is typically taken to be an operator corresponding to a 
conserved current (e.g., particle number or the energy-momentum tensor). While Eq. ([S]) is the Green's function that 
is most closely related to experimental observables, in particular transport properties, it is more convenient in what 
follows to work with the Feynman (time ordered) function, 

iGpiq) = J d4xe*«-(T[O(a;)O(0)]), (6) 

which in the OPE limit takes the form 

The functions Cq(z), with z — /2muj are analytic functions of the variable z, except for single-particle poles at z 1 
and branch cuts at values of z corresponding to multi-particle states. As pointed out in Q, for the theory defined by 
Eq. ([!]), Cq(z) are exactly calculable functions for all z in the case where Oa is in the one or two body sector. If one 
is interested in the asymptotic behavior for — > oo and q — > such that z — )■ 0, the functions Ca{z) can be expanded 
in a power series in z. 

The condensates (Oq) are in general dependent on temperature T, chemical potential /i and the scattering lengtI0 



^ Exactly at the critical point, where the scattering length diverges, the expansion is in fact convergent. 

^ In addition, (Oa) may depend upon other parameters characterizing the ground state, for instance the frequency ujq of a harmonic 
trapping potential. 
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. By dimensional analysis, the functional dependence is of the form 

(0„)=/.^°/2.g„(^^,(aVM)~'), (8) 

so that in the low temperature, a — > cxd limit we can take (Oq) ^ /i'^°/^(7c[(0), and thus the w — > oo limit of a two-point 
correlator in the many-body state becomes an expansion in (non- analytic) powers of /i/w <C 1. Strictly speaking, 
the OPE derived using Eq. ([T]) incorporates only the leading order (scattering length) two-body interaction induced 
by the effective range expansion of the 2 — >■ 2 scattering amplitude. Thus observables calculated using Eq. ([1]) are 
accurate up to corrections parametrized by powers of rovzrriaj, where tq is a length scale characterizing the range of 
the two-body potential (e.g., the van der Waals interaction length scale ivDw)- This means that our results for the 
OPE are valid in the range of energies 

A*<w<-^. (9) 
In this window, the OPE implies the asymptotic (w — > oo, g — )■ 0) expansion 

•=H<,)~^Ev.(9''"(3£:)"'. W 

where ba^i are numerical coefficients proportional to (Oq), and Pa,i are integer. 



B. Sum Rules 



The asymptotic behavior of the Feynman Green's function implied by the OPE, together with analyticity /causality 
properties of the retarded two-point function, can be used to obtain model independent sum rules that must be 
satisfied by the dynamical response functions of the system. It is these analytic properties that will allow us to get 
useful information about the small frequency behavior of the system from the OPE. 

The starting point for this is the Kramers-Kronig dispersion relation satisfied by Gii{q). Because Gii{uj,q) is an 
analytic function in the upper-half complex lo plane Im(cLi) > 0, it follows that 



ReG«(..)=Pr / — 



duj'^ Im Gr{uj') 



Im Gij(w) = -2t^Pr / ' (11) 

These expressions are strictly valid if the function Gr{uj) decays sufficiently rapidly for uj oo. Otherwise, further 
subtraction terms are needed. These introduce additional free parameters and therefore diminish the predictive power 
of the dispersion relations. 

Assuming naively that Gr{lu) decays rapidly enough, it is possible to expand both sides of the dispersion relation 
to obtain, 

Re Gr{oj -^^) = -J2 r^:^ / ^''"Im Gr{u;'). (12) 

It follows from the spectral representation of the two-point functions that, for real to, Re Gr{uj) = Re Gp{ll!). It is 
therefore legitimate to replace the LHS of this equation with Eq. ([7|). If the asymptotic expansion of Re Gi?(cj oo) 
on the real axis consisted of analytic powers of 1/w^, it would then be possible to use Eqs. (|7I12I) to obtain simple 
expressions for the moments 



/•oo T 2 

/„ = / — w2"Im GRiuj'l 
Jo I" 



(13) 



in terms of the condensates {Oa)- Indeed, this is a standard textbook procedure for deriving sum rules. While this 
procedure does lead to several useful moment sum rules (the /-sum rule being among them), it fails but for the first 
few moments /„. The problem is that the asymptotic decay of correlators in any quantum field theory, as predicted 
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by the OPE, is at best power law, and the dispersion integrals in Eq. (jlip are not convergent. 

In order to sidestep this problem, we follow instead the procedure originally used in poj to derive similar sum rules 
in QCD. Consider the contour integral 

I=<f^'-P4, (14) 

taken over a contour consisting of the real axis plus a semi-circular arc of radius R oo through the upper-half plane 
containing the point z = iuj (Ima; = 0,cli > 0). Since Re Gr^—z) = Re Gii{z), and Im Gji{^z) = — Im Gr{z), for 
real z, the residue theorem gives 

doj'^ Im GrIuj' 



-Gr{iuj)^Gf{iu:), (15) 



where we have used the fact that Gr{uj) = Gf{^) for complex uj, which follows by comparing the spectral represen- 
tations. The spectral representation for Gf{uj) also implies that Gpiii^) = Gi?(— w^) is a real function for real lo, and 
thus by successive differentiation, we obtain 

duj'^ Im Gr{uj') 1 f d 



- (..'2+c.2)«+i -n! V rf-V GA-u). (16) 
Combining this with the OPE of Eq. ([7]) , written in the schematic form 

GF{-oj')=Y,a^{q){oj^y'', (17) 

a 

with aa{q) fixed in terms of the condensates (Oq), the sum rules become 



This is a sum rule that relates thermodynamic averages on the RHS (the condensates (Oq)) to the function Im Gu{uj) 
on the LHS. By standard Kubo formulae, Im Gr{uj) controls the dissipative response to perturbations generated by the 
source O and can be measured experimentally (in principle) . Thus the result in Eq. (fT8|) can be either tested against 
experimental results, or used to constrain the parameters of models that make specific predictions for Im Gr{uj). 

For sufficiently low n, the integrals on the LHS of this formula may still fail to converge. However, it is possible to 
take the limit n — > oo, with ujq = uP' jn fixed on both sides of this equation, which yields 



1 (°° du? ^2,2^ ^ , ^ ^ 



r(a) 



K)-"- (19) 



This result should be of phenomenological relevance, provided it is possible to choose ujq large enough so that the 
OPE is valid, yet small enough so that the integral on the LHS is still dominated by the low energy data. 

III. ASYMPTOTICS OF THE DYNAMIC STRUCTURE FACTOR 

We know apply the general results of the previous section to the case of the asymptotic structure factor S{u!,q), 
which is related to the retarded correlator of the densitjd operator 0{x) ~ n{x) = [ip^ 'iIj]{x) , 

Siuj,q) = --lmGR{LO,q), (20) 

TT 

at real frequency w > 0. The asymptotic behavior of this correlator for fermions strictly at the unitarity limit has 
been studied previously in [l3| using the OPE, and by jllj using the methods of ref. [l). 



In what follows, we employ the convention that suppressed spin indices are summed over. E.g., i/jti/i = ^ V'iLV'ai etc. 



6 





FIG. 1. The exact 2 — 2 scattering amplitude. 



Here we will determine the structure of the OPE for 

M{q) = / d*xe''«-^T [n{x)n{0)] 



Ca(g)Oo(0). 



(21) 



Since the OPE holds inside all matrix elements, one may calculate the low-lying Wilson coefficients {q) by evaluating 
both sides of the above expression between states of small particle number, without reference to the interacting many- 
body ground state of the system. Once the coefficients are fixed in this way, the OPE may be evaluated in the 
many-body ground state at finite T and (or for polarized samples with differing chemical potentials for the 

two spin states). 

We will compute the Wilson coefficients in Eq. ([21]) up to dimension A = 5 (in d = 3 spatial dimensions). In 
particular, this excludes operators built out of six powers or more powers of ip, ip^ (by the particle number symmetry 
-0(3;) — e'^"ip{x), operators with odd powers cannot appear). This implies that to fix the relevant coefficients, it is 
sufficient to compute matrix elements of Eq. (|21l) between one-particle and two-particle states. In the 1 — )• 1' and 
2 — > 2' sectors, the theory in Eq. ([Ij is simple enough that the coefficients Ca{q) are calculable exactly. 

In order to compute the relevant matrix elements, we will use dimensional regularization to regulate infrared (IR) 
and ultraviolet (UV) divergences that arise at intermediate stages in the calculations. Our results for the Wilson 
coefficients below will be manifestly IR finite, as well as UV finite after introducing suitably renormalized operators. 

The basic element in the calculation is the exact amplitude for 2 — > 2 scattering, which is given by the sum over 
the graphs in Fig. [T] The bubble expansion in Fig. [T]is a geometric series which sums to 



A'\io,q) = fico.q), 

9 

where {uj,q) is the total energy/momentum of the incoming particles, and 



;T{1 - d/2) [^m{uj - + ie)f^^'^ 



(22) 



(23) 



is the result of computing a one-loop Feynman diagram. In general, this one-loop graph is UV divergent, but in c? > 2 
this divergence is a power divergence that dimensional regularization sets to zero. For example, in d = 3 where the 
one-loop amplitude is linearly divergent, the above result becomes, in the center-of-mass frame. 



47r 

A = 



1 



m 1/a + iy/mto 



(24) 



where the scattering length is 



mg 
4^' 



(25) 



so that in dimensional regularization, the unitary limit corresponds to g — > cxi. 
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FIG. 2. Diagrams contributing to the matrix element of the operator M{q) between one-particle states. The symbol (g) denotes 
an insertion of the operator tp^ili. 



A. Matching in the one-particle sector 



The matrix element of the RHS of Eq. (1211) between one-particle states with momentum p^^ — {p^,p) is given by 
the graphs in Fig. [2] The result, which is exact to all orders in the coupling is given by 

{p'\Miq)\p) = \ + — ^— . (26) 

For reasons that will become clear below, we keep the one-particle states off-shell, meaning that we do not assume 
p'^ — Ep = p^/2m. We have suppressed the spin labels on the external states in this calculation, but note that the 
matrix element is zero unless initial and final spins are equal. 



In order to match onto the OPE prediction, we take the limit g — >■ c» in this expression and expand the non- 
relativistic propagators. Away from the one-particle thresholds q*^ = lu ^ ^E^, the expansion is reproduced by the 
1 — > 1 matrix elements of the operator n(0) = at dimension A = 3, the A = 4 operators 

J- = -2^^^^ V-, (27) 

and the A = 5 operators 



5,9jn(0),a,J,(0) + a,J,(0), 

+ (28) 



The OPE in the one-body sector thus reads 



M{q)\^^y ^ Cn{q)n{0) + Cf\q)d,n{Q) + Cf ^(g) J,(0) + C\f {q)d,djn{Q) 
+ c'ffisD [9.J,(0) -t- d,jm] + (9)^^^ 9.9, V 

^2 

+C(4)(g)V'^ {^^t + ^i: + C(''\q)i^^[-^'^t+^^, (29) 
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with, 



Cniq) 



{uj - E^Yioj + E^Y m' 

AiujE^ 
' {u - E^Yiu + E^Y'^' 



Am? 



1 



E^Y 



1 



' 2m 



{u:~E^Y {^ + E,-Y 
1 1 



{lu - E,-Y (w + Eg^Y 



[uj-E^Y' 
i 

'{uJ + E^Y' 



Note that by the Heisenberg equations of motion for the field operator ip{x), 



2m J 



i, = -25(4vj)(^t^^) 



(30) 

(31) 
(32) 
(33) 
(34) 
(35) 
(36) 
(37) 

(38) 



However, because we are matching to the OPE coefficients using off-sheh states, it is not permissible to apply equations 
of motion inside matrix elements at this point in the calculation. The reason for matching to the OPE using off-shell 
states will become clear below, when we compute the terms in the OPE that are necessary to reproduce the 2—5-2' 
matrix elements of M{q). 



B. Matching in the two-particle sector 



The 2 — 2 matrix element of the LHS of Eq. ((2T|) consists of the graphs in Figs. I3I4I In order to match onto 
operators of dimension A < 5 in the OPE, it is sufficient to compute the (off-shell) two-particle matrix elements 
taking the limit of vanishing external energy and momentum, in which case the calculations simplify considerably. 







FIG. 4. All loop graph topologies contributing to the the 2 — 2 matrix element of the operator M{q). The shaded blobs denote 
exact vertices, given in Fig. [1] The symbol denotes an insertion of the operator ip'^ip. 



The graphs in Fig. [31 with insertions of V'^V' purely on external lines, are given by 



Fig.Ha) =4z^(q) 
Fig. m,b) = 4^(0) 



and 



Fig. He) = 4a(0) 



(39) 
(40) 

(41) 



In this last equation fc = (fc'', fc) — >■ is the energy-momentum of any of the external states. 

Note that the graph in Fig. |3{c) are singular when the external states are taken on-shell, fc° — )■ E^. To regularize 
this singularity, we have kept all external states off-shell in the computation of matrix elements of the bifocal operator 
M{q). For consistency, we have done this even in the 1 — 1' matching calculation of the previous section. This is 
the reason why it is not valid to eliminate the operators in the last line of Eq. (|29p in favor of the quartic operator 
■0|'0|'0TV'i- The 1/0 on-sheU singularities in the 2 — > 2' matrix element of M{q) are reproduced by similar singularities 
on the RHS of the OPE (see Eq. ([5^ below), and cancel out in the computation of the Wilson coefficients. 

Loop corrections to the matrix element are computed by first performing the integration over loop energy by 
residues. The result is, 



Fig. Ha) 



AiA{'d)A{q) 



h{q) + {q^ -q), 



(42) 



for graphs with a single operator insertion on an internal line. The graphs with two internal line operator insertions 
are 



Fig.|l&) = -^(0)2/2 (q) 
Fig.|lc) = ^^(0)2/3 (q) 



(q 
(q 



-q'h 
-q), 



Fig.Hd) = -iA{OfA{q)lUq) + {q^ -q)- 



(43) 

(44) 
(45) 
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In these equations we have introduced 



and, 



These integrals can be performed by a combination of partial fraction decomposition and Feynman parameter tech- 
niques, with the results (with z = E^/{2uj + ie)) 

Ii{q) = -^^(1 - zf/^-^F.ih 2 - d/2; d/2, z/{l - z)), (48) 
Lu + te 

Hq) - - ,'^^i^:°l (l - zY'^-SF,{2, 3 - d/2; d/2, z/(l - z)), (49) 
(cj + ze)^ 

^ , , _ 2/1 (g) , 4m^ r^(d/2 - l)r(2 - d/2) |gf 

where f{u!,q) is the function defined in Eq. ([221) • 

The sum over the above diagrams gives the exact matrix element of J d'^xe^'^'^T[n{x)n{0)] between two particle 
states. In order to match this to the OPE, it is necessary to compare this matrix element to the two-particle matrix 
elements of the RHS of Eq. ([29]). We find that that the matrix element in the q ^ oo limit is reproduced by the OPE, 
provided that it includes the quartic operator 

O4 = V^V'jV't^i, (51) 

with a Wilson coefficient given by 

C^.{q){2'\0,\2) = {2'\M{q)\2) - {2'\J2Ca'{q)Oa,\2), (52) 

a' 

where in the second line the sum is over the one-body operators appearing in Eqs. pop -p7 p . In dimensional regular- 
ization, for off-shell external states, (Oeom = (^idt + "0) 

(2'K0)|2)^4z^(0).^^^, (53) 

k 

(2' I Oeom (0)1 2) = -4^(0), (54) 
(2'|O4(0)|2) = -AHO)/g^ (55) 

while all other bilinear operators have vanishing two-particle matrix elements at zero external momentunj^- We 
therefore obtain. 



Aq) [h{q)--- lihiq) + h{q)) + {q ^ -q) 



(56) 



(57) 

LO — Eg' U! + Eg' 



Note that the potential 1/0 singularities from the on-shell limit of Eq. (HIT) have cancelled against similar divergences 
in the matrix element of Eq. ([55)) leaving behind a finite remainder. We have also checked, using the method 
of regions (l9| . that all IR divergences in the integrals /i,2.3('?) cancel, for arbitrary dimension d, against similar IR 
divergences in the two-particle matrix elements of the bilinear operators in Eqs ([53]) - ([55)) . Thus the Wilson coefficients 
are manifestly IR finite as expected on general grounds. 



More precisely, the matrix elements have power IR divergences for d > 2 which are set to zero in dimensional regularization. 
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Of the seven operators with A < 5, only ri(0), Oij = didj i', and the two operators on the last line of 
Eq. ((29)) . and O4 have, by translation and rotational invariance, non-zero values in the many body ground state. 
Using rotational/translation invariance, (Oij) has VEV proportional to the kinetic energy 



while, by use of the equations of motion 

(^^ (idt dt +1^)) = -25(04). 

Thus for g — > C!0, the OPE predicts 

{M{q ^ 00)) ^ C„(g)(n) - Icff (^IVV^P) + [c^^q) ~ ^gC^'Hq) ~ 25C(^)(9)] 
This is more conveniently written in an operator basis consisting of the Hamiltonian density 

and the contact operator = g'^'>p\'>pl'>P^'>Pi- The result is 



{Miq)) = Cr.iq){n)+Cn{q){n)+CM{Oc), 
where Cn{q) is given in Eq. ((30|) . while 

4i 



Cniq) 



1 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 



and, 



CM = ^ 



Aiq)[hiq)--- — 

g ui — 



- l^h{q) - ^h[q) + {q ^ -q) 



9 d 



1 



(w - Sg-)3 {lo + E^f 



2i 
9 



E? 



(64) 



As a consistency check of these results, we note that the Ward identities implied by current conservation, dtn+V-J = 0, 
impose the constraint M{u!,q — >■ 0) = 0. It is straightforward to verify that the Wilson coefficients in Eq. (15^ all 
vanish at the point g = as required by the Ward identity. We also note that the Ward identity requires matrix 
elements of the operator n{x), and thus M{q) to be renormahzation group (RG) invariant. This property is reflected 
in our result, involving UV finite (in d > 2) Wilson coefficients multiplying the RG invariant operators n, H, and Oj^. 

In fact all three operators n, Oc appearing in the OPE of T[n{x)n{Q)\ have a direct thermodynamic interpretation 
when evaluated in the many-body ground state, in terms of the particle number density N/V, the energy density 
E/V and Tan's contact parameter C, respectively. The latter can be expressed as the thermodynamic derivative. 



C = {0,)=kBTg^dglnZ\^^^ 
which in c? = 3 can be written in terms of the free energy density — F/V 



^ _ 47r dT 



(65) 



(66) 



T,N 



^ RG invariance of Oc is derived in ref . [3| . 



12 



C. Three spatial dimensions 



We now apply these results to the case d = 3 of primary physical interest. In the limit w — oo, z = E^j^uj ^ 1, 
the Wilson coefficients can be expanded as 



cj \rau: 



(67) 
(68) 



and, keeping terms up to order l/.g^ in the 5 — > oo limit 

2 



4? 

Cc{q) = - 



16 2 

-{f{u:) + f{-u)) + --g- 



■9 



(69) 



where f{uj) = —m^—m{oj + ie)/47r. Since for real uj > 0, Ini Gi?(a;) 
length/large frequency limit 



9 V/(^) /(-^) 

Ini Gji{uj), we obtain in the long wave- 



^ \ 2 

q 



- + -(aV^) 



C, 



(70) 



which is valid up to corrections to the OPE from operators of dimension A > 5, as well as corrections from the finite 
range of the two-body potential (neglected in Eq. ([T|). The a 00 limit of our result agrees with the calculation 
of [l3], which uses an equivalent method to ours, but disagrees with that of [ll|. The scaling S{uj,q) ^ w^^/^ is 
consistent with the power counting argument of 21] for the case of a contact interaction. 

Because S{uj,q) falls off as a power for — > 00, the naive sum rules that follow from expanding the dispersion 
integrals of Eq. (ITU) are not applicable. Nevertheless, calculating KeGF{u!,q) for real u > and keeping only the 
leading order term in the OPE, proportional to one recovers the usual /-sum rule 
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Jo 2m V 



(71) 



Perhaps more useful for comparison with experimental data is the Borel sum rule of Eq. ()19p , which requires evaluation 
of the Wilson coefficients at imaginary uj — iujQ, (wq > 0). From Eq. 



Cc{iiOo) 
and therefore, from Eq. ([19] 
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which is the main result of this section. 



IV. POLARIZATION OBSERVABLES 

The methods of the previous section can also be used to study the asymptotic behavior of the polarized density 
correlator 



(74) 
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In particular, the 1 — )■ 1 matrix elements are now identically zero, so that the OPE in this case does not generate 
bilinear operators. In the 2 — >■ 2' sector, the matrix elements are given by graphs with the same topologies as in Figs. 
3,4. Matching with zero-momentum, spin singlet external states, the non-zero graphs are 



Fig.Ha) =^^(g) 



1 2 



+ {q^ -q), 



Fig.EJfe) = 2U(0) 
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Fig.Ha) 

Fig. He) 
Fig.Hd) 

Neglecting operators with A > 5, the OPE for A^^{q) is then 

Mniq) ^ C:-^(g)a(0), 

where now 

2 1 " 
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c^Hq) = 



Aiq) h{q)^ 
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-h{q) + {q^-q) 
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(79) 
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(81) 



Note that due to SU{2) spin symmetry the spin densities n^, ni are separately conserved. The Ward identities 
generated by these conservation laws imply, as in the unpolarized case, that {Mfi{uj,q — >■ 0)) 0. This is reflected 
in our result for the OPE, since it follows from Eq. ([8T|) that C^^^uj, g — >■ 0) = 0. 

As before, the asymptotic behavior of the correlator is determined by condensates that encode the thermodynamic 
properties of the many-body state. Unlike the unpolarized results of the previous section, the asymptotics is now 
dominated by the two-body operator Oc and therefore the Tan contact parameter C. 

For real w > 0, we find, in the q — > limit 



= -^I- G]t(-^^<f) - -{^ (mc.)-3/2C + 0(a-2) 

and thus S-\-i{uj, q) ^ lo^^/"^ in the cj — > oo limit. In addition, we obtain the sum rule for g — > 
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Unlike the sum rule for the density structure function in Eq. (j73p . this sum rule is dominated by contributions from 
intermediate two-particle states and is therefore governed by the behavior of the contact C near unitarity. 



V. BOSONS NEAR UNITARITY 



The results discussed above apply also to bosonic systems near unitarity. A representative example is IIe-4, with 
a ^ 20iydw- Thus to a good approximation, it is described by the Lagrangian, 

^=^'G^*+£)^~i'^''- ^^^^ 

Then the calculations of the previous sections go through virtually unchanged. The two-particle scattering amplitude 
is now 

A-HLO,q) = -j-^f{uj,q) (85) 
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and thus the scattcrmg length m d = 3 is a = mA/STr. In the one-body sector, the OPE of M{q), defined as in Eq. (PT|) 
with n{x) — Lp^Lp{x) is identical to Eq. (|29l) with the replacements tp ^ (p. The 2 — >■ 2' matrix element of M{q) is 
given by the same graphs as in Figs. [31 21 These can be obtained from the results of Eqs. (j5^ - (P5|) by inserting a 
symmetry factor of 1/2 for every loop, and replacing the amplitude by Eq. ([55]) . Using the equations of motion, 



(86) 



the OPE can again be written in the basis of RG invariant operators n{x), the Hamiltonian Hix), and the quartic 

A2 



The result isEI 

{M{q ^ c^)) ^ Cn{q){n) + C,{q){0,), 
where C„ (q) and C-u {q) are given by Eq. ((30)) and Eq. (|63)) respectively, while Cc {q) now becomes 

2 



(87) 



(88) 



lA{q){h{q)-i 



1 



' X d 



1 



1 



- \l2{q) - ^hiq) + {q ^ -q) 



2i 
J 



1 



oj — E,r ljJ + EfT 



{lo - E^f {uj + Eijf 
lin the 

S{uj,q) 



From these results we obtain the asymptotic behavior of the dynamic structure factor 
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as well as the sum rule 
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(89) 



(90) 
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Note that Eq. (|90l) can be compared directly to inelastic neutron scattering data. In particular, the w"^/^ tail in 
S{uj,q) is known to provide a good fit to data [2l|. We plan to apply the sum rule Eq. (PTj) against the He-4 data in 
future work, where in principle one could fit for the value of the contact parameter. 



VI. CONCLUSIONS AND FURTHER DIRECTIONS 



In this paper we have used the OPE to develop sum rules for the dynamic structure function. The sum rule relates a 
weighted integral of 5(0;, q) to a. sum over expectation values of local operators. There are two sources for corrections 
to the rule. First of all the sum over operators is truncated in an expansion in powers of h/ivq, with ^ the chemical 
potential and ujq is a free parameter. There are also small corrections due to the finite range of the two-body potential, 
as well as contributions from operators beyond those with A < 5 considered in this paper. 

It would be interesting to use our sum rules to extract the contact parameter C in various systems where the finite 
scattering length universality corrections become non- negligible. Of course, the accuracy of the extraction is limited 
by the size of finite range effects, which in principle can be included in our sum rules, at least perturbatively. 

In QCD, sum rules similar to the ones presented here have been used to infer properties of low lying excitations. 
This is done by assuming that at low energies the spectral density is dominated by a low lying resonance in a particular 
channel. The continuum contribution has little effect on the values of the low energy parameters. Applying similar 
ideas to cold atoms systems seems like a natural step. 



Note that we have dropped the contribution of the dimension-five operator H. The reason is that for bosons, there is a three-body 
operator of the same scaling dimensions which would need to be included as well [22| . 
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